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ABSTRACT: We study the wash-out of a cosmological baryon asymmetry produced via leptogen- 
esis by subsequent interactions. Therefore we focus on a scenario in which a lepton asymmetry is 
established in the out-of-equilibrium decays of the next-to-lightest right-handed neutrino. We apply 
the full classical Boltzmann equations without the assumption of kinetic equilibrium and including 
all quantum statistical factors to calculate the wash-out of the lepton asymmetry by interactions 
of the lightest right-handed state. We include scattering processes with top quarks in our analysis. 
This is of particular interest since the wash-out is enhanced by scatterings and the use of mode 
equations with quantum statistical distribution functions. In this way we provide a restriction on 
the parameter space for this scenario. 
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1. Introduction and conventions 

Relying on the see-saw mechanism [1-5] to create small masses for the Standard Model (SM) 
neutrinos, leptogenesis [6] provides an attractive explanation of the baryon asymmetry of the uni- 
verse. The main ingredients of leptogenesis are the following: the out-of-equilibrium decays of the 
heavy right-handed neutrinos, that were added in the see-saw mechanism, into leptons and Higgs 
particles violate CP, from which a lepton asymmetry can be generated. This lepton asymmetry 
is then partially transformed into a baryon asymmetry by anomalous processes of the SM called 
sphalerons [7]. In this way the three Sakharov conditions are fulfilled and leptogenesis turns out as 
a consequence of the see-saw mechanism. 

In the minimal (type-I) see-saw model, gauge singlet fermions with Majorana masses Mj, 
eigenvalues of the complex symmetric 3x3 Majorana mass matrix Mm, couple to the massless 
lepton doublet and to the Higgs doublet of the SM through Yukawa couplings given by the 3x3 
matrix \ v . These Yukawa couplings then generate a Dirac mass term = Xv, where v = 
174GeV is the vacuum expectation value of the Higgs field, relating the heavy singlets to SM 
neutrinos upon spontaneous symmetry breaking of the electroweak gauge symmetry. In the see- 
saw mechanism the eigenvalues of the Majorana mass matrix are much larger than the Dirac mass 
matrix entries, Mm, > m^, and the mass matrix of the left-handed states is then suppressed by 
the high-energy scale M, 

m Vij = -v 2 X ik M M l k \ kj . (1.1) 
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Participating in Yukawa couplings, the heavy neutrinos are unstable and the total decay width of a 
right-handed neutrino generation TVj at tree-level is given by 

r A =Y l (r ia (N i ^$ + l a )+T ia (N i ^* + l a )) = (At A) Mi, (1.2) 

a 

where the sum is taken over the single decay rates into lepton flavor a. The CP asymmetry in lepton 
flavor a is then defined as 

-l ia lex) Tja 

This asymmetry is due to the interference of the tree-level amplitude and the one-loop vertex and 
self-energy contributions. The indices i and a denote the generation of the decaying right-handed 
neutrino and the flavor of the produced lepton, respectively. Accounting for different lepton gen- 
erations, the CP asymmetry is a diagonal matrix in flavor space. It has been shown in [8-12] that 
the flavor structure of the final state leptons has sizable implications on leptogenesis predictions. 
To account for the lepton flavor structure one introduces the projectors 

AP° 

P ia ^\(k\l a )\ 2 = P? a + ^, (1.4) 

APO 

P ia = \(l'i\l a )\ 2 = P^-^. (1.5) 

The wash-out rate in inverse decays is then reduced by the projector Pf a = (Pj a + Pi a ) /2, since 
the Higgs will make inverse decays on flavor eigenstates \l a ) instead of the linear superposition 
\li). The second effect, an additional CP-violating contribution, stemming from a different flavor 
composition between \U) and CP |7j), can be expressed by the projector difference APj a = Pj a — 
Pi a . The sum over flavor of this difference results in J2 a AP® a = 0. The decay rates into single 
lepton flavor a are now given by Tj Q = Pi a Ti and Tj a = Pi a Ti, and the flavored CP asymmetries 
can be written in terms of the unflavored asymmetries 

A P 

£ia = Pia Si + — g— ■ (1.6) 

The second term gives the contribution due to a potentially different flavor composition between 
\li) and |IA. Summing these CP asymmetries over the flavor yields the total CP asymmetry used 
in the unflavored regime E; L = J2 a Ei a . In this paper we will mainly be interested in the effects 
of a full spectral distribution of right-handed neutrinos on leptogenesis and not in the effects of 
the lepton flavor structure. Therefore we will perform all calculations in the so-called alignment 
case, that is realized when the Ni -decays are just into one flavor a, resulting in Pi a = Pi a = 1, 
Pi0^ a = Pip^a = and ei a = E{. In this case the kinetic equations, c.f. Section 3, are identical to 
the kinetic equations used in the unflavored regime that is realized at T > 10 12 GeV. By restricting 
to the alignment case we can compare our results to the ones obtained in previous calculations 
using unflavored Boltzmann equations [21-23]. However, we will keep all flavor indices in our 
notation, since all results are valid when considering independently the single flavors a. We will 
comment on the effects of the lepton flavor structure in more detail later in Section 2. 
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Summing over all lepton flavors, = ^ £j Q , the CP asymmetry can be written as [13-15] 



£,; — 



- E 



Im 



(AtA) 



/MIA /Mf 



(1.7) 



where the vertex (_/y) and the self-energy (/s) contributions are 

1 + x 



fv(x) 



1 - (1 + x) log 



, and f s (x) 



(1.8) 



The out-of-equilibrium dynamic that is necessary for successful leptogenesis is provided by the 
expansion of the universe. One usually compares the total decay rate of the right-handed neutrino 
state given in Eq. (1.2) to the expansion rate at temperatures T ~ Mi, 



H(T = Mi) = V47r3 ff */45 (M,/M P1 ) , 



introducing the decay parameters 



K - TDi 
1 ~ H(Mi 



mi 



(1.9) 



(1.10) 



where M Pi = 1.221 x 10 19 GeV is the Planck mass, and g* = 106.75 corresponds to the number 
of relativistic degrees of freedom in the SM at temperatures higher than the electroweak scale. One 
introduces at this point two dimensionful variables to connect the decay parameters to the neutrino 
mass scale via the effective neutrino mass [16] fhi and the equilibrium neutrino mass m*: 



rru 



and m* 



16712-^/77* v 



Mi ' " ' 3^/5 Mpi 
The decay parameters for each flavor a are again given via the projectors 

K = P° K 

— tot - fv «- 



1.08 x 10" 3 eV. 



(1.11) 



(1.12) 



The decay parameters control whether the right-handed neutrino decays in equilibrium (Ki a > 1) 
or out of equilibrium (Ki a < 1) and are key quantities for the dynamics of leptogenesis. 
In general, the baryon asymmetry produced by leptogenesis can be written as [17] 



3 a sp h 



£ i* K L = d^e ia n{ a ~ 0.96 x 10 2 ^ e { k: 



/ 



(1.13) 



where in the last step the sum over the flavor a has been performed as it is possible in the alignment 
case or in the unflavored regime at T > 10 12 GeV. The final efficiency factors k{, that parametrize 
the amount of asymmetry that survives the competing production and wash-out processes are direct 
results of solving the relevant Boltzmann equations for leptogenesis [18, 19]. In the limit of vanish- 
ing wash-out and a thermal initial abundance for the right-handed neutrinos, the efficiency factors 
have a final value k{ = 1. The factor / = 2387/86 accounts for the dilution of the baryon asym- 
metry due to photon production from the onset of leptogenesis till recombination and the quantity 
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a S ph = 28/79 is the conversion factor of the B — L asymmetry into a baryon asymmetry by the 
sphaleron processes. 

The most simple realization of thermal leptogenesis is the N\ -dominated scenario, where a 
lepton asymmetry is exclusively generated in the decays of the lightest right-handed state N\. 
Qualitatively, this scenario can be realized assuming M\ <C Mi <C 10 14 GeV <C M3. Then the 
heaviest right-handed state N3 decouples [20] and this typically implies a hierarchy in the CP 
asymmetries, i.e., \es\ <C \s2\ <C |ei|, as a consequence of light particles in the internal lines 
of the self-energy and vertex corrections to the tree-level decay diagram. Indeed, in the limit of 
massless particles in the internal lines, the corresponding CP asymmetry vanishes, cf. Eq. (1.7) . 
Furthermore, due to N\ interactions following the decays of the heavier right-handed states ^2,3, 
a substantial part of the produced lepton asymmetry will presumably be washed-out, especially in 
the regime of strong wash-out. 

In this paper we will go beyond the simplest scenarios and consider leptogenesis in a set-up 
where the asymmetry is generated in the decays of the next-to-lightest state N2, that is called the 
^-dominated scenario [21]. The viability of ^-dominated leptogenesis is restricted by the wash- 
out of produced asymmetry induced by interactions of the lightest state N\ that tend to destroy 
any previously created lepton asymmetry. We will address this issue of wash-out in ^-dominated 
leptogenesis within the simple alignment case, that is equivalent to the unflavored regime, by means 
of complete Boltzmann equations at the mode level including scattering with top quarks [22] and 
extend on the previous studies of [23] where only decays and inverse decays were considered. We 
show how the additional wash-out due to quantum statistical distribution functions and scattering 
processes restricts the valid parameter space of ^-dominated leptogenesis. 

In the next section we show how specific scenarios of leptogenesis can be realized by specify- 
ing the relevant see-saw parameters and discuss qualitatively the N\- and ^-dominated scenarios. 
In Section 3 we present Boltzmann equations for leptogenesis including scattering at the mode 
level and apply these equations to study ^-dominated leptogenesis quantitatively in Section 4. In 
the Appendices A, B, and C we list explicitly the integral expressions used to solve the Boltzmann 
equations. 



2. The orthogonal matrix and different scenarios of leptogenesis 

To understand how specific scenarios of leptogenesis can be realized, we recast the light neutrino 
mass matrix, Eq. (1.1), 

1 T 

m v = -m D ——m D , 

M M 

with mo = Aw. Here, it is always possible to choose a basis in which the heavy neutrino mass 
matrix is diagonal, Dm = diag (Mi, M2, M3). Using an unitary matrix U, one can simultaneously 
diagonalize the light neutrino mass matrix 

D m = -U^m u U*, (2.1) 

where D m = diag (mi, 7712,7/13). If one does not account for the running of neutrino parame- 
ters from the electroweak scale to the see-saw scale [24, 25], the matrix U corresponds to the 
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Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix. In the basis where the charged leptons mass 
matrix is diagonalb and with the help of an orthogonal matrix 17, the Dirac neutrino mass matrix 
can be written in the so-called Casas-lbarra parametrization [26], 



mo 



(2.2) 



The Dirac neutrino mass matrix is fully described by 18 parameters: the mixing matrix U contains 
six parameters (three mixing angles and three phases), the diagonal matrices D m and Dm contain 
three neutrino masses each, and the orthogonal matrix Q is described by six real (three complex) 
parameters. It can be written as a product of three rotational matrices [21,27] 



^ (^21, W31, LO32) — Rl2 (W21) Rl3 (W31) R23 (^32) 



(2.3) 



with 



R12 — 
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and 



R23 = 



Rl3 — 



±1 



-W31 
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"^32 







W32 



^32/ 



(2.5) 



(2.6) 



In general, one can state that Eq. (2.2) is divided into two parts: (i) a measurable low-energy part, 
containing the PMNS matrix U and the light neutrino masses D m , and (ii) a high-scale part, con- 
sisting of the orthogonal matrix and the heavy neutrino masses Dm, which is not accessible by 
current experiments. 

The Ni -dominated scenario can now be realized assuming the heavy neutrino mass matrix to 
be hierarchical, i.e. Mi <C M 2 <C M 3 , together with a specific choice of the CP asymmetries 
implemented by a special form of the 17 matrix [21,27]: 

• For Q, = R13, the CP asymmetry in N2 decays vanishes, i.e., 62 = 0, while £\ is maximal. 

• For $7 = R12, the CP asymmetry e\ is suppressed compared to its maximal value, and 
\e2\ oc (M1/M2) |ei| is negligible within a strong mass hierarchy x . 

In [17,28,29] it has been shown that the N\ -dominated scenario proves to be independent of the 
initial conditions in the regime of strong wash-out. However, in the weak wash-out regime the final 



'if Mi ~ Mi, both CP asymmetries should be taken into account. 
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asymmetry production depends sensibly on the initial conditions of the N\ abundance. Further- 
more, thermal leptogenesis sets a lower limit of Mi > 10 9 GeV [30] on the mass of the lightest 
right-handed neutrino in order to explain successfully the observed value of the baryon asymmetry 
of the universe. This bound is consequently translated into a lower bound on the reheating temper- 
ature after inflation, Trh ^ 1-5 x 10 9 GeV [31]. These lower bounds not only have an issue with 
the cosmological abundance of gravitinos [32-35] within supersymmetric theories, but also cause 
some specific problems in grand unified theories based on flavor models. Some of these models as- 
sume a grand unified symmetry between up-quarks and neutrinos. The neutrino Yukawa couplings 
are then connected with the up-quark Yukawa matrices leading to right-handed neutrino masses 
that are proportional to the square of the up-quark masses [36, 37]. Typical values for the mass 
of the lightest right-handed state fall in the range 10 6 — 10 7 GeV, see e.g., [38, 39], and masses 
Mi > 10 9 GeV need a specific choice of parameters [40,41]. This makes thermal leptogenesis 
difficult to reconcile with this class of models. 

In order to circumvent these issues, the ^-dominated scenario was proposed in [21]. Indeed, 
for Q, = i?23 one can have a maximal CP asymmetry S2 coming from N2 decays while e\ van- 
ishes. This means that N\ is totally decoupled from the heavier states while in the N2 decays the 
heavy third state, N%, runs in the internal lines of the one-loop diagrams. Together with a mass 
hierarchy in the heavy neutrinos, ^-dominated leptogenesis can be realized if the wash-out from 
Ni interactions does not deplete the produced asymmetry. With the above choice of $7, the effec- 
tive neutrino mass of the lightest right-handed state is fixed to m.\ = m\. Thus, for hierarchical 
light neutrino masses the N\ interactions can be forced to be in the weak wash-out regime where 
rhi < m* ^ 10~ 3 eV. It is worth noticing that, if the N\ interactions are in the weak wash- 
out regime, i.e., rh\ < m*, then the N2 interactions are constrained to be in the strong wash-out 
regime, i.e., rfi2 > m*. This is due to the orthogonality of the Q matrix [21]. Therefore, N2- 
dominated leptogenesis is independent of the initial conditions on N2. In order to achieve a large 
enough CP asymmetry, the lower bound on the mass Mi can be directly translated into a bound 
on the mass M2, whereas the lower bound on the reheating temperature remains unchanged. The 
lower bound on Mi is then evaded and Mi can be arbitrarily light compatibly with the see-saw 
condition M S> mo- However, allowing for small complex rotations R12 and i?i3, both of these 
bounds become increasingly more stringent leading to a point beyond which the ^-dominated 
scenario is not viable anymore. 

Note on flavor 

Considering flavor effects beyond the alignment case may substantially change the parameter 
ranges in which the scenarios discussed above are valid. When the flavor structure of the lepton 
asymmetry and the wash-out is tracked, it is possible to generate a large enough lepton asymmetry 
in N2 decays even if Ni interactions are effective. It has been shown in [41] that the Ni -induced 
wash-out of an asymmetry, established in the lepton state I2 via N2 decays 2 , can be sizably ef- 
fected by the projector P\ a = |(/i|^ Q )| 2 that leads to a reduction of the wash-out strength, i.e., 
K\ a < K\. Another possible modification is that the coherence of the lepton state I2 is destroyed 
by potentially fast occurring N\ inverse decays. This leads to a statistically mixed state consisting 

2 The states \h) = l/y/ (AtA) i; A a »|2a) are participating in interactions with the heavy states Ni [45]. 
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Table 1: Scenarios in the decays/inverse decays pictures (D 1 and D4) and scenarios including scattering with 
the top quark (SI and S2) and their associated assumptions. Cases Dl and SI correspond to the conventional 
integrated Boltzmann approach, while Cases D4 and S2 involve solving the full set of Boltzmann equations 
at the mode level. The naming of the different scenarios corresponds to [22]. 



Case 


kinetic equilibrium 


quantum statistics 


Dl 


Yes 


No 


D4 


No 


Yes 


SI 


Yes 


No 


S2 


No 


Yes 



of l\ and the state orthogonal to l\. As a result, parts of the asymmetry, stored in the state orthogonal 
to li, are protected from Ni wash-out [42-44]. This is only possible when the Ni states decay in 
the two-flavor regime at temperatures 10 9 GeV <T< 10 12 GeV, where only the r Yukawa inter- 
actions are in equilibrium, or in the unflavored regime at temperatures T > 10 12 GeV. In the fully 
flavored regime, however, the full flavor basis is resolved and there does not exist any direction 
in flavor space that is protected against Ni wash-out [45]. The study of [46] considered the case 
U = R12 (W21) R13 (^31) together with W32 = for values Mi > 10 9 GeV, lying in the two-flavor 
regime. It was found that even for small values u 2 i < 0.1 the contribution to the asymmetry from 
N 2 decays can be naturally the dominant one for 10 11 GeV < M 2 < 10 14 GeV, when the effects 
of flavor on the N\ -induced wash-out are taken into account. 

Furthermore, in Section 1 we did not account for the effect of flavor on the CP asymmetry. 
In Eq. (1.3) we wrote the CP asymmetry as a matrix in flavor space. Indeed, the hierarchy of CP 
asymmetries, following a mass hierarchy of particles in the internal lines, does not necessarily hold 
when flavored asymmetries are considered [11]. In contrast to the unflavored scenario, where e 2 is 
suppressed by M\jM 2 compared to E\ when assuming M\ <C M 2 <C 10 14 GeV <C M3, there is 
the possibility of having a non-negligible asymmetry e 2a that, in turn, eventually extends the range 
of the ^-dominated scenario. 



3. Full Boltzmann equations for leptogenesis 

In this paper we study the wash-out of a lepton asymmetry by means of the full Boltzmann 
equations at the mode level including scattering with the top quark in the alignment case where 
Ki a = Ki and K ia ^ = 0. This study does not rely on the assumption of kinetic equilibrium 
for right-handed neutrinos and will not restrict the distribution functions to be classical. In order 
to study the effects of quantum statistical distribution functions and scattering with top quarks we 
specify different scenarios shown in Table 1. The mode equations including scattering processes 
based on the Yukawa coupling with the top quark for right-handed neutrinos and the lepton asym- 
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metry are given by 

H(Mi) df Ni 



H(Mi) d f(i-i) a 



= C D [f Ni ] + 2C s , s lf Ni ]+4Cs,t [f Ni ], 



dzi 



C D 



■fy-I) J + 2C S,s [f(l-l)a\ + 4C S,t 



(3.1) 



(3-2) 



where z.- L = Mi/T. The right-hand sides of Eqs. (3.1) and (3.2) contain one collision integral from 
decays/inverse decays and two collision integrals from scattering processes coming, respectively, 
from scattering in the s-channel and in the t-channel at tree level. Here, one factor of 2 stems from 
contributions from processes involving anti-particles, and another factor of 2 in the t-channel terms 
originates from the u-channel diagram. 

In our analysis we only include tree-level scattering processes that are of O (h 2 X 2 ) [22]. We 
do not consider interactions with gauge bosons, nor include CP violation in 2 — > 2 or 1(2) — > 3 
processes, which are of higher order in the Yukawa couplings. The effect of CP violation from 
these processes was considered in [9,47-49]. The numerical integration of Eqs. (3.1) and (3.2) over 
the right-handed neutrino (lepton) phase space results in the time evolutions of the number densities 
njVj and n^_ ^ . Normalization to the equilibrium photon number density n 7 q = (C(3)/7r 2 )g 7 T 3 , 
with g 7 = 2, gives the comoving abundances (N^_j^ ). The decay collision integrals for the 
right-handed neutrinos and the lepton asymmetry can be evaluated analytically to be [23,50] 



C D 



Cd [In, 



£nvn 



sinh((£jv - y N )/2) 



(3-3) 



where F+ 



+ o 

fl a + fl a > £ 3 



N 



\ 2 \ 2 



f(i-T 

: Ej/T, Vj 



_sinh((£:Ar + yjv)/2) 
(/* + f*i)(f(i-T) a + e <« F « ) " 2e - M (i + /*) 

(3.4) 



Pj|/T, 3 and /$ denotes the distribution function of 
the Higgs particle. Up to subleading corrections, that are of higher order in the Yukawa couplings, 
the collision integrals in the different flavors decouple and the corresponding Boltzmann equations 
can be solved in analogy to the single flavor case. 4 For SM particles we assume Fermi-Dirac distri- 
butions for leptons and a Bose-Einstein distribution for the Higgs particle, since gauge interactions 
lead to very fast equilibration rates [51]. 

The scattering collision integrals are nine-dimensional and cannot be evaluated analytically. 
However, they can be reduced analytically to two-dimensional integrals given by [22] 

3T h 2 t {T)Mimi 



C S,(s,t) [fNi] = ^2 



2 6 7r 3 £ N y N 



x 



N,l,q,t 



d£i 



w{£ N ,l,q,t) 



k(£ 



N,l,t,q 



d£ (t,i) A (^) 7 (£t)' 



(3-5) 



(3-6) 



3 For better clearness we drop the generation index for the heavy neutrinos in En and j/jv and the flavor index for the 
leptons in yi . 

4 Reference [23] includes terms in addition to Eq. (3.4) in order to avoid asymmetry production in thermal equilibrium. 
However, the same analysis also shows that the quantitative difference between this and our approach is negligible. 
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for the right-handed neutrino and 



C S ,(s,t) [/(,_!) J -E^3^2 



3T h\(T)Mim x 



v 2 



ru(£ N ,i, q ,t) rk(£ N ,i,t, q ) , - n , . 

< / ^ / A (3-7) 

J W(tN,l,q,t) J l(t-N,l,q,t > 



for the asymmetry. Here, /?| (T) is the top Yukawa coupling evaluated at the leptogenesis scale 
and the functions ^ and ( P a denote the phase space factors for the right-handed neutrino 
(lepton) scatterings in the s- and i-channel, respectively. The integration limits u, w, k, and I depend 
on the energies of the particles involved in the interactions and give the distinct integration ranges 
for each integrand in the sums of Eqs. (3.5) and (3.7). The integrands consist of the phase space 
factors and the analytical functions ^ , that encode the matrix element and the result of the angle 
integration. Both depend on the energies of the particles involved in the process. For the s-channel 
diagram the sum contains six terms for the right-handed neutrino and the lepton asymmetry and for 
the ^-channel diagram one counts four terms each. The explicit expressions of the single collision 
integrals are listed in the Appendices A and B and the rather lengthy derivation can be found in the 
appendices of [22,52]. 

Inserting Eqs. (3.5) and (3.7) in the Boltzmann equations (3.1) and (3.2) yields the complete 
set of differential equations to be integrated numerically. Once the Boltzmann equations for the 
distribution functions are solved, one can perform the integration over the right-handed neutrino 
(lepton) phase space to obtain the number densities and, in turn, the efficiency factors 

4 er 1 

3 a sp h - 1 " l ) a 

A value of -^(;_j) = 3.14 x 1CT 8 is needed in order to explain the value r/§ MB = (6.225 ± 

0.17) x 10~ 10 , measured in Cosmic Microwave Background (CMB) observations by the WMAP 
satellite [53]. 

With the restriction on classical equilibrium distribution functions of the Maxwell-Boltzmann 
type and assuming kinetic equilibrium holds for the heavy right-handed neutrinos, Eqs. (3.1) and (3.2) 
can be integrated over the right-handed neutrino (lepton) phase space to recast the Boltzmann equa- 
tions in an integrated form [17], 



dN Ni 
dz 



= - (A + Si) [N Ni - N c ^j , (3.9) 
-Wia N^_- l)a + e ia Di (aTjv, - , (3.10) 



where Nj = rij/n^. The quantities Di {Si) account for the production of right-handed neutrinos 
from inverse decays (scatterings), and the wash-out rate Wi a contains as well contributions from 
both, scatterings (W<f ) and inverse decays W?® = Pf a The explicit expressions for Di, Si, 

Wi a , and N^. can be found in Appendix C. 
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Figure 1: The final efficiency factors Kj a at the production stage for a vanishing initial Ni abundance in 
dependence of the interaction parameter Ki a . Shown are the two cases including scattering processes S 1 
(dotted/blue) and S2 (short dashed/magenta), and two scenarios Dl (solid/red) and D4 (long dashed/green) 
within the decay-inverse decay only framework [22]. 



4. Mode equations in A^ 2 -dominated leptogenesis 

The general effects of the full spectral distributions of right handed neutrinos on the production of a 
lepton asymmetry in thermal leptogenesis have been discussed extensively in [22], where the final 
efficiency factors for the scenarios, presented in Table 1, have been calculated. The results can be 
summarized in Figure 1, where the final efficiency factors as a function of Ki a are shown for the 
integrated approach (Cases Dl and SI) and the complete mode treatment (Cases D4 and S2), both 
including and excluding scattering. For Cases SI and S2 which include scattering, we note that 
their difference is rather large in the weak wash-out regime (Ki a < l),with the integrated approach 
overestimating k{ q by up to a factor ~ 1.5 at Ki a ~ 0.01 compared to solving the complete mode 
equations. But this difference decreases as we increase K{ a . At Ki a > 3 the integrated approach 
underestimates n ia by less than ~ 10%. It is also interesting to note that the relative contribution 
of scattering processes to the final efficiency factor is smaller in the complete mode calculation 
than in the integrated approach. In the weak wash-out regime, including scattering enhances the 
final efficiency factor from decays and inverse decays by up to a factor of ~ 30 in the integrated 
scenario. In the complete mode calculation, however, the enhancement is only a factor of ~ 15. 
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Similarly, in the strong wash-out regime, scattering reduces n{ a by up to 20% in the integrated 
picture, compared to below 10% in the complete treatment. 

In this section, however, we want to discuss the effects the complete set of mode equations, cf. 
Eqs. (3.1) and (3.2), have in the iV2 -dominated scenario. 

Here, the wash-out due to N\ interactions is (i) enhanced using mode-equations with de- 
cays/inverse decays alone and (ii) enhanced including scatterings with the top quark. We will not 



0.001 




Figure 2: Lepton asymmetry generated in N2 decays for e^a = 10 6 (M2/IO 10 GeV) and different values 
of M2 in the Case S2 for a vanishing initial N2 abundance: (i) M2 = 10 13 GeV, solid line, (ii) M2 — 
10 12 GeV, long-dashed line, (iii) M 2 = 10 n GeV, dashed line, (iv) M 2 = 10 10 GeV, point-point-dashed 
line, and (v) M2 — 1.9 x 10 9 GeV, point-dashed line. Within the considered scenario, the shaded region 
indicates where the produced baryon asymmetry exceeds ?y§ MB . The allowed values of jiV/^lin will be 
taken as initial values to calculate the wash-out induced by N\ -interactions. The differences induced when 
considering the Cases Dl, D4 and SI can directly be read-off from Figure 1. 



study ./V2 -dominated leptogenesis with a specific form for the Q matrix but rather consider values 
of K\ a and K2 a going from the weak wash-out regime (Ki a = 0.01) up to the strong wash-out 
regime (Ki a = 100) without the restrictions that are usually imposed when choosing a specific 
matrix constellation. For example one usually has K 2a > -^soi ~ 9 for O23 [21,46]. In doing so 
we can study the effect of the enhanced wash-out due to scattering and the use of mode equations 
on the iVi -induced wash-out and restrict the general available parameter space for ^-dominated 
leptogenesis in terms of and K^ 2 )a- A similar study has been performed in [23] where, 
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however, wash-out effects have been considered in the decay-inverse decay only picture. 

The lepton asymmetry that can be generated in N2 decays is shown in Figure 2 in dependence 
of the decay parameter Ki a for values of M 2 varying between 1.9 x 10 9 GeV and 10 13 GeV. The 
shaded area marks the region where the final asymmetry exceeds the value ??§ MB , measured in the 
observation of the Cosmic Microwave Background. For M 2 < 2x 10 9 GeV the CP asymmetry 
generated in N2 decays is too small to account for the observed value of the matter-antimatter 
asymmetry. 5 We assumed the maximal value of the CP asymmetry e™ x ps 10 -6 (M2/IO 10 GeV) 
to be realized. In specific scenarios, however, the CP asymmetry E2 a will be smaller and the 
corresponding lower bound will be at M 2 > 10 10 GeV [21]. On the other hand, for M 2 > 10 13 GeV 
one has to account for AL = 2 violating scattering processes with the N2 in the s- and t-channel. 
Being an additional contribution to the wash-out, these processes tend to reduce the final amount 
of asymmetry [54]. The two vertical lines in Figure 2 correspond to the solar and atmospheric 
neutrino scales, respectively. Considering the N2 interactions to be in the strong wash-out regime 
in the window preferred by neutrino oscillation data demands M2 > 10 11 GeV to explain the 
observed value of the baryon asymmetry. 

The lepton asymmetry generated in N 2 decays is altered by the subsequent wash-out due to 
interactions of the lightest right-handed neutrino N±. According to the considerations of [23], where 
./V2 -dominated leptogenesis has been addressed by means of mode equations within the decay- 
inverse decay only scenario (Case D4), we choose the following initial conditions at z\ = M±/T to 
calculate the effect of wash-out on an initially produced asymmetry: (i) We take |AT j|| in = 1CT 7 
as initial value of the lepton asymmetry generated in 7V 2 decays, (ii) assume a zero initial N\ 
abundance, and (iii) set e± a ps 0. The third condition can be achieved by supposing a small value 
of the Ni mass. Anyway, small lepton asymmetries stemming from different generations add up 
linearly and an additional asymmetry e\ a would not modify our considerations on the N\ -induced 
wash-out effects. We choose M\ = 10 7 GeV in the numerical implementation in order to fix the 
evolution of the top Yukawa coupling. 

Figure 3 shows on the left panel the time evolution of the normalized N\ number density for 
the Cases Dl, D4, SI, and S2 in dependence of z\ and M\ = 10 7 GeV. On the right panel the time 
evolution of the lepton asymmetry during Ni wash-out is shown for the same scenarios. Thus, in 
addition to the discussion in [23], we include Cases SI and S2 here. In the weak wash-out regime 
(Ki a = 0.1), the asymmetry is only slightly reduced in the Cases D4, SI, and S2 compared to the 
integrated approach in the decay-inverse decay only scenario, Case Dl. The net wash-out of the 
initial asymmetry is less than 10%. However, already at K\ a ~ 1 the strength of the wash-out in the 
different scenarios becomes distinguishable. At z\ ~ 1 wash-out becomes effective and is strongest 
in Case S2 where the complete set of Boltzmann equations including scattering with the top quark 
is considered. However, the net reduction of the initial asymmetry is still less than one order of 
magnitude. The difference in the lepton asymmetry between Case D4 and Case S2, both using 
mode equations, is very small (below 2 %) and cannot be see in Figure 3 . 6 Comparing Case S2 
with Case SI, we see that the influence of the additional wash-out factor f^ 1 , present in the mode 
equation, cf. Eq. (3.4), is larger than the wash-out due to scatterings in the integrated approach that 

5 This bound on M2 is more severe than the value given in Section 2 since for a vanishing initial abundance n{ a < 1. 
6 An error in the numerical calculation for the Case D4, appearing in an earlier publication [52], that leads to larger 
differences between the Cases D4 and S2, has been corrected in the present paper. 
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K la = 0.1 K la = 0.1 




Figure 3: Time evolution of the absolute value of the normalized right-handed neutrino number density 
Nm 1 /N^ l and the lepton asymmetry \Nq_-^ | for three different coupling strengths K\ a in the N 2 - 

dominated scenario for an initial asymmetry | ^V^^. | i n = 1CP 7 , AI\ = 10 7 GeV and a vanishing initial 
Ni abundance. The differences between the Cases D4 and S2 are smaller than 2 %. 

is present in Case SI. In the strong wash-out regime, for K\ a = 5, the initial asymmetry is depleted 
by up to two orders of magnitude, with the strongest wash-out again in Case S2. Considering the 
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Figure 4: Amount of initial asymmetry that has to be generated in N2 decays in order to survive the subse- 
quent wash-out by iVi interactions for the Cases Dl, D4, SI and S2. The mass of the lightest right-handed 
neutrino was set to Mi = 10 7 GeV and the CP asymmetry generated in Ni decays was set e\ a = 0. The 
mass M^ 2 , that is needed to produce enough initial asymmetry in iV2 decays when K2 a ~ 1, is shown on 
the right vertical axis for the Case S2. The asymmetry in the gray shaded region cannot be generated in N2 
decays and the area right of the arrow is excluded due to N\ wash-out. 



momentum integrated scenarios, the reduction in Case SI is two orders of magnitude larger than 
in Case Dl. Concerning the scenarios in which mode equations are used, the contribution of the 
scatterings amounts again less than 2 % as can be seen by comparing Case S2 with Case D4. In 
general, it can be stated that the wash-out is enhanced more by the use of mode equations (due to an 
additional wash-out factor f^) than by the inclusion of scattering processes that become important 
only for values K\ a > 1 in the integrated approach (Case SI). 

Figure 4 shows the amount of initial asymmetry generated in N2 decays that is needed to 
account for the observed value of the baryon asymmetry after wash-out due to N\ interactions 
in dependence of K\ a for the Cases Dl, D4, SI and S2. For the Case S2, the mass of the next- 
to-lightest neutrino M% 2 , that is needed to produce enough initial asymmetry when Kia. ~ 1> is 
shown on the right vertical axis. The green colored regions correspond to the value of the baryon 
asymmetry deduced in Big Bang Nucleosynthesis (BBN) within 95% confidence level [55], and the 
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dashed-lines within this regions represent the value of the baryon asymmetry deduced from CMB 
measurements. When Ni interactions fall into the weak-wash-out regime {K\ a < 1), almost all 
of the initial generated asymmetry survives giving |-/V^| in ps 3 x 1CT 8 , i.e., the same limit as can 
be seen in Figure 2 and the different Cases Dl, D4, SI and S2 cannot be distinguished. However, 
increasing the strength of the N\ interactions, wash-out becomes more and more effective and for 
Kia ^ 1 the differences between the Cases Dl, D4 and S2 become visible. Again, considering the 
treatment with mode equations (Cases D4 and S2), the contributions of the scatterings are below 
2 % and too small to be visible in Figure 4. For K loi = 10 an initial value |AT ji| in ps 1 is needed to 
account for the observed value of the baryon asymmetry. Though, as can be seen in Figure 2, values 
of the asymmetry lying in the shaded area above the horizontal dashed line at | AT^| in = 1.7 x 10~ 4 
cannot be generated in N2 decays. Therefore, for the ./V2 -dominated scenario to be successful, 
the N\ interactions are restricted to K\ a < 5 for the Cases S2 and D4. When considering the 
integrated Case Dl, this bound can be relaxed to K\ a < 7. Considering Case S2, the upper bound 
on the initial asymmetry corresponds to values Mf 2 = 10 13 GeV and K<i a ~ 1 as can be seen on 
the right vertical axis of Figure 4. For values of K2 a lying in the strong wash-out regime, that is 
preferred by neutrino oscillation data, the generated asymmetry is roughly one order of magnitude 
smaller, leading to K\ a < 3 in Case S2 and K\ a < 4 for Case Dl. When choosing conservative 
values, Mi ps 10 11 GeV and k 2oi ~ 1CP 2 , the scenarios are forced to K\ a < 2 for Case S2 and 
Ki a < 3 for Dl. 7 This corresponds to values of K\ a typically needed in 50(10) inspired GUTs 
where flavor effects are included. For the limiting scenario, M2 ~ 2 x 10 9 GeV, the N\ interactions 
strength K\ a has to vanish [56] in all considered cases. 

5. Conclusions 

In this paper we studied leptogenesis in an alternative scenario where the lepton asymmetry is 
created in the decays of the next-to-lightest right-handed neutrino state N 2 . Here, the additional 
wash-out present in the complete set of mode equations leads to a more efficient depletion of the 
lepton asymmetry in interactions of the lightest right-handed neutrino N\ ; these interactions follow 
the asymmetry generation in the decays of the heavier state. In order to account for the observed 
value of the matter-antimatter asymmetry, the possible values of the decay parameters K\ a and 
K2a of the two right-handed states can be restricted. From the maximal amount of asymmetry that 
is achievable in N2 decays, the decay parameter of the lightest right-handed neutrino is forced to 
Ki a ^ 5 when accounting for the spectral distribution of right-handed neutrinos. The effects of 
scattering process on the N\ -induced wash-out are smaller than 2 %. Furthermore, demanding the 
decay parameter K2 a to be in the strong wash-out regime favored by neutrino oscillation data, 
where the asymmetry generation is independent of the initial conditions on N2, sets the more 
stringent limit K\ a < 2. 

Note added 

After the present paper was submitted, a new analysis of leptogenesis in the ^-dominated scenario 
appeared [58], in which the effects of flavor on the CP asymmetries in ^-dominated leptogenesis 
are studied in detail. 

7 Using typical assumptions, the study in [23] found a limit K\ a < 3. 
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Appendices 

A. s-channel collision integrals 
Right-handed neutrino 

Using energy conservation, and Fermi-Dirac statistics for the leptons and quarks, the phase space 
factor reads 

e Si+£t (-\ + f N . +e £ » f N .) 



4 Nt) (fNi,fl a ,ft, fq) (1 + &£l) (1 + eft) ^ £i+£n +e£t y 

The s-channel contribution to Eq. (3.5) is then given by, 

CsM = + cp + cp + c?> + c& + 



(A.l) 



(A.2) 



The single collision integrals must be evaluated numerically and then summed to give Cs, s [fNi 
The integrals cjg' s '"' 6 ^ are as follows: 

1. First integral (with q = q/T): 
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(£n-Vn) (2£i+gjv-j/jv) 
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2. Second integral: 
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3. Third integral: 
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4. Fourth integral: 
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5. Fifth integral: 
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6. Sixth integral: 
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Lepton asymmetry 

For the lepton asymmetry the s-channel phase space element is given by 

(i-i) / \ e £t (1 + ( e £l+£N - 1) /jv ) 

As " [f(i-T) a >f =-f(i-i) a { i + e s t){e e l+ s N+e e t) » (A. 15) 

and the s-channel contribution to Eq. (3.7) can be expressed as 

cs„lfy- n j = eg + eg + eg + eg + eg + eg <a.i 6) 

to be integrated numerically over two remaining degrees of freedom. The explicit integrals in 
Eq. (A. 16) are: 

1 . First integral 



(1) 3r hlMjrhi fVW^f ,„ f £ i+ £ 
s > s 2Wtf v 2 



/ d£ N / d£ t At l) 4 \ (A. 17) 



where 1^ is given by Eq. (A.4). 



(2) 

2. Second integral with Is given by Eq. (A. 6): 

r (2) 3T fegM.m, [°° [ £ > +£ » (2) 

C ^ = Mf^^J V ^ d£N J, d£tK ° ^ ■ (A " 18) 

3. Third integral (7 S (3) given by Eq. (A.8)): 

4. Fourth integral (ll 4 ^ given by Eq. (A. 10)): 
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5. Fifth integral (/| 5) given by Eq. (A. 12)): 
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6. Sixth integral (ii 6) given by Eq. (A. 14)): 
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B. i-channel collision integrals 
Right-handed neutrino 

The t-channel phase space element for right-handed neutrinos is given by: 

e si+e q (_ l + fNi+e s NfN ^ 



(fNi,fq, flcft) 



(1 + e £ i) (1 + e £ i) (e £ i + e £ i +£ i 
and the t-channel contribution of Eq. (3.5) can equivalently be written as 

CsAfN i \ = c ( H + cfl + c ( il + c ( il 

where the constituent integrals are as follows: 
1. First integral (with k = k/T): 
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2. Second integral: 
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4. Fourth integral: 
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In the integrals (B.3) and (B.4) we have introduced ah = m^/Mi as an infrared cut-off for the 
t-channel diagram, where m<j> is the mass of the Higgs boson which presumably receives contribu- 
tions from interactions with the thermal bath, i.e., its value does not correspond to that potentially 
measured at the LHC. The value of m$> can in principle be deduced from a thermal field theoretic 
treatment of leptogenesis, and the analysis of [28] found m$(T) ~ 0.4 T. However some open 
questions still remain and hence in the present work we prefer to adopt a value of ah = 1CT 5 , used 
first by Luty in [57]. 



Lepton asymmetry 

The t-channel phase space element for the lepton asymmetry is given as 

\( l ~ J ) a ( f f f f \_ f e£g ( e£l tllMz^lM 
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As in Section B we cut-off the integrand in and by adding in the lower integration 
limit of k ( modifying the limits of S q and S n accordingly) . 

The t-channel contribution of the collision integral (3.7) can then equivalently be written as 

Cs,Af(i-l)\ = C { sl + Cfl + eg + eg, (B.8) 
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2. Second integral: 
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4. Fourth integral: 
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C. Reaction rates in the integrated picture 

The decay interaction term Di in the Boltzmann equations (3.9) and (340) is given by 



Di = Zi Ki 



/Ml 

\E N 



(C4) 



The scattering rate Si consists of two terms, Si = 2S S + 4S t , coming respectively from 
scattering in the s-channel and in the i-channel. One factor of 2 stems from contributions from 
processes involving anti-particles, and another factor of 2 in the t-channel term originates from the 
u-channel diagram. 

It is convenient to rewrite the s- and t-channel scattering rates S^ s ^. in terms of the functions 
/(«,*)«» defined as 
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such that 
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K Si 
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and the total scattering rate is given by 



where 
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with mj and m* defined in Eqs. (1.11). 
The functions Xs,t(%) are defined as 
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where again = m^/Mi is the infrared cut-off for the i-channel diagram. 
The total wash-out rate is given by 



Wt = W IDi 



1 + ^( 2 ^ 5 - +4 ^ 



where Wj £> s quantifies the strength of the wash-out due to inverse decays 

1 N N- 
W ID . = -Di -4rs, 

and the equilibrium abundances are A^ eq = 3/4, and Nff (z) = § z 2 Ki (zi). 
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